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I. INTRODUCTION 

Everyone knows the popular child puzzle called "15". It consists of a 4 x 4 grid of 
numbered sliding square tiles with one free slot (vacancy). Initially the tiles are jumbled. In 
the course of the game, not only the empty slot moves around the grid, but the (marked) 
tiles also move, and can take different relative positions for a fixed position of the vacancy. 
The goal of the game is to unjumble the tiles by successive moves, each one consisting in 
sliding a tile into the empty slot. It is easy to see, and this is true for a grid of arbitrary 
size, that the vacancy can freely reach an arbitrary position. 

To complicate the game, take now rectangular tiles -dimers or dominoes- of size 2 x 1 or 
1x2 instead of square ones, and consider tiling a large grid with those, leaving, as before, 
an empty slot of unit surface, i.e. a monomer (so the number of sites in the grid must be 
odd). The possible movements of the vacancy are now conditioned by the orientation of the 
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neighboring dimers, since a dimer can slide only if it is oriented towards the vacancy. The 
main question, raised and studied in is then: what are the positions that the vacancy 
can possibly reach_^under these moves ? 

Bouttier et al [1] have analyzed this question by looking at various quantities related to 
the statistics of the domain accessible to the vacancy (unlike in the game "15", the dimers 
are not marked). Building on the well-known Temperley correspondence between fully- 
packed dimer configurations (i.e. with no vacancy) and spanning trees, they showed that 
the correspondence can be extended to dimer configurations with a single vacancy, which are 
then put in bijection with spanning webs. Generically, a spanning web consists of a central 
tree, growing from the site where the vacancy is located, surrounded by a number of nested 
loops to which further branches are attached. They proved that the domain accessible to 
the vacancy is exactly given by the tree component. Given a uniform distribution on the 
dimer configurations for a fixed position of the vacancy, the previous correspondence induces 
a distribution on the spanning webs with statistical weights depending on their number of 
loops, and implicitly defines a non-trivial distribution on the clusters of sites accessible to 
the vacancy. 

Putting the vacancy at the central site of a square, odd-by-odd grid, Bouttier et al 
performed numerical calculations on finite grids and then extrapolated their findings to very 
large or infinite lattices. One of the main results of their analysis is that the delocalization 
probability Vl-, defined as the probability that the cluster of sites accessible to the vacancy 
covers the whole grid, decays like for large L, the linear size of the grid. From this 

result, the probability p{s) that the domain accessible to the vacancy has size s, in infinite 
volume, was estimated to decay like s~^/® for large s. At the other end of the scale, the 
values of Pl(s) for small values of s were computed exactly for finite grids, and used to 
assess the asymptotic values in infinite volume. The numerical values reported in [1] for 
p(l), the probability that the vacancy be strictly jammed, and p(2), the probability that 
the vacancy can make exactly one move, are p(l) = 0.107864 376 269 049 51198(1) and 
p(2) = 0.055 905 353 801942(1). Moreover, using Plouffe's inverter (2| applied to the first 
ten digits of 1/ a/p(1), the authors of [1] were able to conjecture the exact value of p(l), 
namely p{\) = 57/4 — 10^/2. The exact value of p{2) remained unknown. 

Our aim of this paper is to revisit these questions, and to provide an analytical derivation 
of some of the results mentioned above, mainly the values of p{l) and p(2), as well as other 
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quantities, which are relevant to the diffusion properties of the vacancy. In particular we 
obtain the exact value of p{2), 

p{2) = ^{72 817^2 - 102 977), (1.1) 

which agrees with all significant digits given in |l|. 

The plan of the article is as follows. In Section [III we review the bijection between dimer 
configurations with the fixed vacancy and spanning webs, as formulated and developed in 
|l|, and discuss the spanning web enumeration technique. In Section UTTl we discuss the 
situation when the spanning webs can include a maximal, finite number of loops, in the 
infinite volume limit, and relate it to the finite-size scaling of the delocalization probability 
found in jl|. In Section [IV] we analyze the probability that the vacancy take at least one 
step in a particular direction, while Section |V] contains the results announced above for p{l) 
and p{2). 



II. THE MODEL AND THE SPANNING WEB REPRESENTATION 

Consider a square grid C of odd size. Let the vertices be colored in black and white like 
a chessboard, so that neighboring vertices have different colors. For definiteness, we assume 
that the corner vertices of the grid, which all have the same color, are white. With this 
coloring, there is one more white vertex than black ones so that the grid can be fully covered 
by dimers with a single vacancy on a white vertex. Fig{T^ shows a possible configuration of 
dimers with the vacancy (represented as ®). 

According to [1], a dimer lying next to the vacancy can move (by one site) and cover the 
vacancy if it is oriented towards it. After the sliding of the dimer, the vacancy has moved by 
two sites, but remains on the white sublattice. In the example of Fig{T^, the upper and the 
right neighboring dimers of the vacancy can be moved; equivalently, the vacancy can move 
two sites upwards or two sites rightwards. 

The coordinates of each white vertex are either both even or both odd, so there are two 
types of white vertices, even and odd. The sublattices of even and odd white vertices are 
called respectively even and odd as well, and denoted £even and C^dd (-^even is slightly smaller 
than £odd)- If the vacancy or a dimer are located on an even (odd) vertex, then we will call 
them respectively an even (odd) vacancy and an even (odd) dimer. 
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(a) (b) (c) 

FIG. 1: (a) Example of dimer configuration with one odd vacancy (marked by 0) on a 11 x 11 
grid. The other two panels show the corresponding spanning webs defined (b) on the odd sublattice 
Coddj a-iid, (c) on the even sublattice -Ceven- 

It is easy to see that the parity of the vacancy does not change during its motions. Assume 
that it is odd, like in the example of Fig{T^. To set up the correspondence between dimer 
and spanning web configurations, we first remove the black and even vertices and keep the 
odd dimers only. We then replace each odd dimer by an arrow of length 2, starting from 
the odd vertex and directed along the dimer. The arrow configuration obtained from the 
example in FigHK is shown in FigHb. The set of all possible arrow configurations obtained in 
this way from dimer configurations with a fixed odd vacancy forms the set of spanning webs 
defined on the odd sublattice. As illustrated in Fig{T}3, a spanning web consists of a central 
tree component, rooted at the vacancy location and oriented towards it, surrounded by a 
number of nested directed loops with branches attached to the loops. The sites contained in 
the central tree are precisely the sites which are accessible to the vacancy. Indeed, looking 
back at the original dimer configuration, one sees that the central tree is enclosed by a loop 
of even dimers which acts ge for the vacancy. If the central tree covers the whole 

of the odd sublattice, no loop is allowed and the spanning web is a spanning tree. This is 
automatically the case when the vacancy is at a boundary site. 

Conversely, if we have a fixed spanning web on the odd sublattice, we can reconstruct the 
dimer configuration up to the orientation of the dimer loops on the even sublattice. Indeed 
the set of even dimers also form a spanning web on the even sublattice, which is dual to the 
one on the odd sublattice (see FigHt). The central piece of an even spanning web is a loop, 



(a) (b) (c) 

FIG. 2: (a) Example of dimer configuration with even vacancy, (b) Odd spanning web obtained 
from (a), (c) Even spanning web obtained from (a). 

whose interior is the domain accessible to the odd vacancy. The even spanning web has the 
same number of loops as its dual odd web, but its tree components are oriented towards the 
exterior of the lattice. An odd spanning web fixes the shape of the even spanning web but 
not the orientation of its loops (if any). Since there are two possible orientations per loop, 
2"- dimer configurations correspond to each odd (or even) spanning web with n loops. So 
we can enumerate all dimer configurations on the original lattice by counting the spanning 
webs on the odd (or even) sublattice with a weight 2"'. 

The same construction can be made when the vacancy is even (see Figl2]), but in this 
case the situation is slightly asymmetric with respect to the choice of even or odd spanning 
webs. Indeed it is not difficult to see that the odd spanning webs contain one more loop 
than the even spanning webs. This affects the formulae with a factor 1/2 or 2 depending on 
the choice one makes, namely Zdimer = l^oddweb = 2Zevenweb- 

Which spanning webs and vacancy, even or odd, we use is a matter of convenience. In 
the odd vacancy is located at the center of a (4L + 1) x (4L + 1) grid, and the authors 
chose to work with the odd spanning webs. 

We finish this section by recalling how the enumeration can be carried out. We assume 
that the vacancy is odd, and we choose the spanning webs on the even sublattice £even, 
taken to be finite for the moment. As remarked above, the outmost piece of such a web is a 
(disconnected) tree, rooted at the exterior of £even- From Kirchhoff's theorem, the number 



6 



of spanning webs for which the outer tree fully covers £even, i-e. spanning trees, is equal to 

Zt,ee = detA°P, (2.1) 
where A°p is the discrete Laplacian on £even with open boundary conditions, 

{4 ifz=j, 
-1 if i,j are nearest neighbours, (2-2) 
otherwise. 

In order to count not only the spanning trees but also the spanning webs, we have to allow 
for loops circling around the vacancy (which lies on Codd)- This can be done by choosing 
a defect line connecting the vacancy to the boundary, for instance the dotted line shown 
in FiglH:;. Then we define a frustrated Laplacian matrix A°p, equal to A°p except that 
(A°P)jj = +1 for all those nearest neighbour sites connected by a bond (of £even) which 
crosses the defect line. Because any loop enclosing the vacancy necessarily contains an odd 
number of frustrated edges, it can then be shown that the determinant of A°p counts 
the number of spanning webs with the correct weight, namely 2" if the web contains n loops, 



^web 



detA°P. (2.3) 



If {ik,ik) are the pairs of neighbouring sites separated by the defect line, for k = 
1,2, . . . , N, where N depends on the size of the grid, the position of the vacancy and the 
defect line that has been chosen, the matrix A°p can be seen as a perturbation of A°p, 
A°P = A°P + 5, with 

Bi^j^ = Bj^^i^ = 2, Bij = elsewhere. (2.4) 

This allows to relate the number of spanning webs to the number of spanning trees, 

1^ = det(/ + G°P5), (2.5) 

where is the Green matrix, (^"p = (A°p)^^. This number is the inverse of the localization 
probability introduced in ij]. Let us note that the previous formula requires that .^tree = 
det A be non-zero, and so is not defined in the case of an even vacancy if the odd spanning 
webs are used, since A would be subjected to the closed boundary conditions, implying 
det A'^i = 0. 
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III. FINITE DEFECT LINES 



Consider first the situation of a finite defect line, of length A^. The system is first confined 
in a finite grid; the infinite volume limit is then taken, while keeping finite but large. This 
situation is geometrically depicted in FigOl 
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FIG. 3: Typical loop configurations in presence of a defect line of length N. The vacancy (on the 
left) is odd, while the sites ik,jk are even. 

The perturbed Laplacian corresponding to the finite defect line is A°p = A°p + B^^^ with 
the non-zero part of B^^^ given by 
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(3.1) 



0/ 



The defect matrix B^^^ is a 2N x 2N matrix, with indices ranging over zi, . . . , iN,ji, 
We want to compute the ratio 

det A°P 



— 



det(/ + G°P5(°)), 



,JN- 



(3.2) 



det A°P 

which is the fraction of arrow configurations on the even sublattice with at most loops, 
each containing an odd number of frustrated bonds, compared to the number of spanning 
trees. In terms of dimers, we can think of two odd vacancies placed at the two ends of the 
defect line, with the loops enclosing one of the two vacancies. 
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If the even sublattice is a finite grid of size L, the determinants in the numerator and de- 
nominator are x L^, but because the matrix B^^^ has rank 2N, the rightmost determinant 
in f l3.2p is only 2N x 2N. When L goes to infinity, the order of this determinant remains 
equal to 2N, and the inverse Laplacian 6*°^ may be replaced by the inverse Laplacian G on 
the infinite lattice. Then the simple block form of -B^'^^ makes it possible to factorize the 
corresponding determinant into two A^-by-A^ determinants of Toeplitz matrices: 



vP = det 



s 


S' 


S' 


s 



S-S' 


S' 





S + S' 



det{S - S')Ndet{S + S')n, (3.3) 



with the matrices S, S' given by 

Sk,£ = 5fc/ + 2Gj^j^ = + 2Gjj.^i^, (3.4) 

for k,i = 1,2, . . . , N. The translational invariance of the Green function, Gf^^f^' = G^_fi ^ = 

Gp^q with {p, q) = r — f', shows that S and S' are Toeplitz matrices, and the reflection 
symmetries imply 

Sk,e = Sk,e + '^Gk-e,i , and S"^ ^ = 2Gk-e,o ■ (3.6) 
From the explicit formula, 

G„ = r ^ (3.7) 
' 47r ^(2 -cos 1^)2-1 



we find the generating functions 



+ 00 



n=— oo 

+ 00 



- E (S» - S>"- = ^/^^^l^ = yf^. (3.9) 

\'6 — 2cosv9 /+(V5) 



The determinants of S* ± 5" can be computed by using Widom's theorem, a generalization 
of Szego's theorem j^. Let f{<f) = X] '^ne'"'^ be a function on the unit circle of the form 

n 

f{if) = {2-2cos<frg{<f), «>-^' (3-10) 

where g{^p) is a smooth univalent function, nowhere vanishing nor divergent. Then the 
asymptotic value of the Toeplitz determinant formed with the Fourier coefficients a„ of /, 
is given by 

D^_i(/) = det{ae.k)i<k,e<N ^ E[g- a] N"" e^^^^^^)", N » 1, (3.11) 
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where (log (7)0 is the zeroth Fourier coefficient of log (7, and E[g; a] is a constant whose explicit 



value can be found in j4|. 

Applying this theorem to /+ and /_ yields 

{\ogg+)o = — / log(/+(<^)rf<^ = log(l + y2), (log^7_)o = -log(l + v^), (3.12) 

Jo 

and consequently, 

= det{S ± S')n c:^ E± m [l + V2j , (3.13) 
where E± = E[g±; =l=|]. Finally, we obtain the ratio (13. 2p . for large N: 

pSJ^ ~ A^i (3.14) 

We should notice that E+ is actually divergent, with a divergence proportional to Go.o- 
This is clearly due to the fact that /+ has a non-integrable singularity {a^ = —h), which man- 
ifests itself in J'j^ by the constant E[g+; a+] being singular when — > — | Physically, 
this is not surprising either, since there are much more (infinitely more in thermodynamic 
limit) spanning webs than spanning trees. However the divergence does not depend on N, 
which means that we get the same "infinity" for = 1 as for any other value of A^. This 
suggests the better definition of V^^^ as the relative fraction of spanning webs with at most 
A^ loops with respect, not to the spanning trees, but to the spanning webs with at most 
one loop. The corresponding ratio, after a proper regularization (see |5|] for instance), would 
then be well-defined. In any case, we will ignore the divergence, since only ratios of P^"* are 
relevant for what follows. 

The quantity V^^^ as such does not quite compare with the delocalization probability 
studied by Bouttier et a/ in [l|. We may however, cut the plane by inserting a vertical 
boundary passing through the sites ii and ji, and consider the left half-plane only. In this 
way, the arrow configurations forming spanning webs with loops around the single vacancy 
are generated, and we recover a situation similar the one studied in jl|, except that we 
work on a half-plane rather than on a large square. Which boundary conditions, open or 
closed, are imposed on the vertical boundary make no difference. The above formulas (13.31) . 
(13. 4p and (13.51) remain valid provided one uses the Green function appropriate to the chosen 
boundary conditions. The resulting S, S' matrices in (13. 6p are no longer Toeplitz but pick 
a Hankel piece, depending on k + £. We have nevertheless computed the corresponding 
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determinants numerically. Our numerical results for this case show that the power 1/2 in 
Eq. (13.141) is changed to 1/4, which agrees with the exponent found in [l|. 

It is instructive also to compare the asymptotics (13.141) with that obtained by Fisher and 
Stephenson [6] for the monomer-monomer correlations. Considering a pair of monomers 
located at two sites separated by an odd number of bonds of the original lattice, they ob- 
tained the asymptotic correlation proportional to N~^^'^, with the distance between the 
monomers. In our case, the vacancies belong to the same sublattice and are therefore sepa- 
rated by an even number of bonds. The difference between the two asymptotics demonstrates 
the crucial role of the loop statistics which is quite different in these cases. 



IV. MOVING ONE STEP 



As a preparation, we start by computing the probability Vmob that the vacancy can move 
at least one step in a specific direction, in the infinite volume limit. Due to the rotational 
symmetry, Vmob can be interpreted as an efficient mobility of the vacancy in any direction. 
For convenience, we use in this section the spanning webs defined on the odd sublattice, so 
that the sites ik,jk defining the defect line are on the same sublattice as the odd vacancy, 
located at the site in+i, see FigJH 

With respect to the discussion in Section [Tll the difference implied by the use of the odd 
sublattice is that the Laplacian of reference is now A^,, namely the Laplacian with closed 
boundary condition and a root at iN+i- For actual calculations, one can take (A^)^^ = 
(A'^^)jj + e Si^i^^-^ ^j,iN+i which case the number of spanning trees growing from the root 
is equal to lim -A^. As in the previous section, loops are included by using the perturbed 
Laplacian Ag, = A^ + B'^^\ The asymmetry of the descriptions afforded by the even and odd 
vacancies, mentioned in section II, implies that the number of dimer configurations with one 
vacancy is equal Zdimer = hni - det A^ for an odd vacancy or Zdimer = h det A^^ for an even 
vacancy, with A'^^ = A'^' + B^^\ On a finite grid, these two numbers are slightly different, 
but become equal in the thermodynamic limit. 

Suppose that the vacancy can take one step upwards, so that the dimer covering the 
odd site Jn+i is oriented toward in+i, like on FigHJo. In terms of arrows, there is an arrow 
pointing from Jn+i to ^a^+i- The number of arrow configurations with this arrow from Jn+i 
to iN+i is equal to the number of arrow configurations for which Jn+i is also a root (in 



11 



jN+l JN 

O 

e 



iN+\ 



6 
iff 



Jn-1 





6 



(a) 



J2 

o 



l2 



Jl 

Q 



6 
h 



Jn+\ 



FIG. 4: Figure (a) describes the pertubation needed to force the presence of a dimer pointing to 
the vacancy, as shown in (b). 



addition to iiv+i)- Indeed the two roots generate separate branches attached to them, but 
these branches can be joined to a single tree, if we place a fixed arrow (dimer) from Jn+i 
to ijv+i. In this way, we generate all spanning webs on the odd sublattice with an odd 
vacancy at tN+i and an arrow on the bond [ja^+i, ■^at+i]. Therefore, by enumerating these 
spanning webs and dividing the result by the total number of spanning webs without the 
arrow constraint, we obtain the probability Pmob that the vacancy can make a step upwards. 

Turning the site j'^v+i into a root requires to further perturb A|g, by the diagonal term 
^^i,jN+i ^j,jN+i- The resulting Laplacian can be written as A'^^ + B^^\ where the non-zero 
entries of B^^^ are 



B. 



*JV + l!*JV + l 



B 



(4.1) 
(4.2) 



On a finite grid, the probability that the vacancy can make a step upwards is equal to 
the ratio 



lim ^det(A^i + S«) 

£— >CX3 ^ 

lim -detAsa 



(4.3) 



In the thermodynamic limit, the denominator can be replaced by |det A'^^ and A*^' by the 
Laplacian on the plane A, so that the previous ratio goes to 



lim 4 det(A^i + B^^^) lim 4 det(A + B^^^) lim ^ det(J + GB^^^] 



(1) 



det(A'=i + S(o)) 



det(A + 5(o)) 



det(/ + GBW) 



2^- 



V 



{0)-' 



N 



(4.4) 

of which the limit N oo readily yields Vmoh- 

The nonzero part of the defect matrix B^^^ consists of two blocks of size + 1. As in 
the previous case, the corresponding determinant "D^"* can be factorized into two Toeplitz 
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matrices of size + 1, 

p« = -1 detiS -S'- detiS + S' - 5^''+'^)n+u (4.5) 

for the same matrices Sk/ and S'f. ^ as in the previous section, except that their indices run 
from 1 to + 1, and where the matrices 5'-''^ are defined as 6^1 = 6k^s^e,s for any s = 1, 2, . . . . 
The expansion of these determinants by the last element (A^ + 1, + 1) gives 

det(S ±S'- 5^^+^%+! = J-^+i - J^. (4.6) 

Therefore, using (13.13^ . (13.14^ and (14. 4p . we obtain 

2P^^^ 

P^ob= lim ^ = V2-1. (4.7) 

V. LOCALIZATION PROBABILITIES OF THE VACANCY 

Dimer configurations with a vacancy on the odd sublattice are in correspondence with 
spanning webs on the even sublattice. As recalled above, the central part of an even spanning 
web consists of a loop surrounding a certain number of odd sites of the original lattice, which 
exactly form the domain accessible to the vacancy. Following pj], let Pl{s) be the probability 
that this domain has size s when the grid is finite (with size proportional to L), and let p{s) 
be its infinite volume limit value. Hereafter, we compute p(l) and p(2), the probabilities 
that the vacancy can take respectively zero or one move, so that p(l) is the probability that 
the vacancy be fully jammed. We start with p{l). 

If the vacancy is fully jammed, it has to be surrounded by four dimers, arranged as in 
FiglSb or a similar pattern in which the dimer covering points west. In the two cases, 
the arrows drawn on the even sublattice are constrained to form a minimal loop around the 
vacancy, so that the arrows originating from the four sites iN,iN+i, Jn, Jn+i niust form a 
cycle. 

Let us now turn the sites inAn+i^ Jn^ Jn+i into roots, like in the previous section. On the 
rest of the even sublattice, the arrows are all unconstrained, and can possible point towards 
iN,iN+i, Jn or Jn+i- If in addition, we insert a defect line of length as before, the arrows 
may form loops passing through the bonds (zi, ji), {12,32), ■ ■ ■ (^TV-i, jAf-i)- These changes 
lead to the perturbed Laplacian = A°p + B^'^\ where the non-zero entries of the defect 
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FIG. 5: In (a) is a pictorial description of the perturbation needed to account for the presence of 
a minimal loop surrounding the vacancy, as shown in (b). 



matrix are given by 



i-k,3k 3k,^k 



IJV,«JV ijv+l,«JV+l 



?{2) 



JN,JN 



1,2,. ..,N, 

JN+1,JN+1 ' 



(5.1) 
(5.2) 



for e very large. 

The determinant lim ^ det counts the number of even spanning webs with at most 
— 1 loops in presence of the four roots. Every such spanning web gives rise to two dimer 
configurations on the whole lattice with nine vacancies, namely the original vacancy and 
the eight surrounding sites. In turn these eight sites can be covered by four dimers in two 
ways. Thus the number of dimer coverings in which four dimers loop around the vacancy 
(and with the number of even loops bounded by A^) is equal to 4 x lim ^ det A^, and their 
fraction in the set of similar dimer coverings but with no restriction around the vacancy is 



4 X 



lim ^det(A°P + 5(2)) 



4 X 



lim ^det(/ + G'°P5(2)) 



(5.3) 



det(A°P + 5(o)) det(/ + G°P5(o)) 

In the thermodynamic limit, (^"p can be replaced by the Green matrix on the plane, and 
the previous ratio becomes 

1(2) 



-r,y^, lim if det (/ + G5 (2)) 



V 



(0) 
N 



= 4 X 



det(J + 6-^(0)) 



(5.4) 



Finally the limit N ^ oo gives the jamming probability p{l), 

p(l) 

where P^"* is known from fl3.14p . 



4 lim ^ 



(0)' 
N 



(5.5) 
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(2) 

The rank 27V + 2 determinant defining V)^^ can be factorized as before, 

vP = ^ det(^ - 5' - - 5(^+^))jv+i X det(5 + S' - ^^^^ - 5(^+^))Ar+i, (5.6) 

and yields the following form for the jamming probability p{l), 

1 . det(S-y-5(^)- 5(^+1)^+1 det(S + 5'-5W-5(^+i)W , , 
^(') = 4i^oo detiS-S% ^ d^^{STS% • ^'-'^ 

Expanding the determinants and using the values of J'^ = det(5' ± S')n given in fl3.13p 
recasts the ratios as 



det(^±^'-5(^) _ J^^, Mi^N;N){S ± S') 



det(S±5')iv det(^±5')^ J^"" Jt 



N 



to leading order in A^. The notation Mi(^N;N)iS ± S')n+i stands for the principal minor 
det(S'ij ± Slj)i<ij<N+i,i,jj^N- 

Minors of Toeplitz matrices have been very recently studied by Bump and Diaconis, who 
have obtained general formulae f?! (see also [8|]). The expressions useful for the case at hand, 
and for the evaluation of p{2) below, are briefly recalled in the Appendix. Using them, we 
obtain the asymptotic value 

M-';-^)'^ = 1 + (log/,),.(log/J. ^ 4 - 2V/2. (5.9) 

N^oo det{S ± S')n 

The two ratios (15. 8p are then both equal to 4v^ — 5 for large A^, and yield the result quoted 
in Q, 

1 57 

p(l) = -(AV2 - 5f = - - IOV2. (5.10) 

The similar calculation can be carried out for p{2). If the vacancy can take one and only 
one step, the dimers around it must be organized as in FigjUjD. In terms of arrows drawn on 
the even sublattice, they all correspond to fix a rectangular length 6 loop running around 
the vacancy. The shape of the loop can be rotated by 90 degrees, and each loop may have 
two orientations, making eight different situations. They however all contribute the same 
number. 

The arguments used for p{l) apply to this case as well, leading to a perturbation of the 
Laplacian involving six roots and a defect line of length N. The corresponding defect matrix 
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FIG. 6: Geometric setting for the calculation of p(2), corresponding to a 6-dimer loop around the 
vacancy. 



B^^^ therefore reads 
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(3) 
'ikJk 
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1,2, 



B'^'. =B''" . =5^^^ . =B\ 

^N,^N *JV+l,*iV+l *iV+2,*Ar+2 JN,JN 



?{3) 



JiV+lJJV+1 JJV+2JJV+2 



e. 



(5.11) 
(5.12) 



Reasoning as before, the determinant lim det(A°P + B^^^) counts the number of even 
spanning webs with at most — 1 loops in presence of six roots. Again every such spanning 
web gives rise to two dimer configurations on the whole lattice with fifteen vacancies, namely 
all the sites shown in Figl6]D. If one excepts the vacancy itself, the other fourteen sites can 
be covered by seven dimers in eight different ways, among which only two contain a 6-dimer 
loop. So the number of dimer coverings in which six dimers make a loop around the vacancy 
(and with the number of even loops bounded by A^) is equal to 4 x lim det(A°P + B^^^). 
Their fraction in the set of similar dimer coverings but with no restriction around the vacancy 
is, in the thermodynamic limit, equal to 

lim idet(/ + G5(3)) 
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4^^ 
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det{I + GBio)) 



(5.13) 



Taking the limit N ^ oo and multiplying the result by 4 to account for the four rotations 
of the rectangle yields the jamming probability p{2), 



p{2) 



p(3) 

16 lim ^ 



(0)' 
N 



(5.14) 



The factorized form of T>^m reads 



N 



V 



(3) 



1 



_ _ ^ det (S - - - 5(^+1) - 5(^+2) ) ^+2 X det (5 + - - 5(^+1) - 5(^+2) ) _ 
64 



(5.15) 
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Straightforward algebra reduces the ratios to previous results and two new types of minors, 
det(5±y-5(^) -5(^+1) - 5(^+2))^^^ ri2v^- 16^ 



det(^±^')iv+i I 4v^-8 

'^i{N;N)iS ± S')n+2 , . ^ M^N,N+l;N,N+l)iS±S')N+2 , .^s ^r;ic\ 

- detiS±S%,, + detiS±S% ^ ^ ^ 

to leading order in A^, and where Mi(7v,Ar+i;iv,Ar+i)(5' ± S')n+2 is the determinant det(S'jj ± 
^ij)'i-<i,j<N+2,i,jj^N,N+i- The evaluation of these two minors is again a simple apphcation of 
the Bump-Diaconis formulae j?]. Surprisingly we find that the four minors are equal by 
pairs, 

hm ^^(^'^)'''^ ^ '^'-'^+2 _ j.^ Mi(Ar_Ar+i;Ar^Ar+i)(S' =F 5")7V+2 _ I ^ — 46a/2, (^IJ) 

N^oo det{S±S')N+i ~ N^c>o det(^T^')iv ~ I 73 _ ' 



Inserting these results in the above equations, we obtain the exact value for p{2), 

p(2) = i-(72 817^2 - 102 977), (5.18) 
as it was quoted in the Introduction. 
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APPENDIX: MINORS OF TOEPLITZ MATRICES 

We collect here the various formulae used in Section |V] to compute specific minors of 
Toeplitz matrices. These are but particular cases of more general expressions proved by 
Bump and Diaconis in [7]. In fact the asymptotic values of exactly the same determinants 
have been obtained by Tracy and Widom 81] at about the same time. However their general 
answer takes a different form and appears to be less convenient for concrete calculations. 

Generically these formulae evaluate the determinants of Toeplitz matrices with certain 
rows and columns removed or shifted, and which are, as a consequence, no longer Toeplitz. 
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For f{(p) = a-n^^^'^ a function on the unit circle, let Tjv-i(/) = {a^-k)i<k/<N be the 

n 

usual rank Toeplitz matrix associated to the symbol /, and Dj^_i{f) = detT/v-i(/) its 
determinant. 

Now let A = (Ai, A2, A3, . . .) be a partition of m, that is, a decreasing sequence of 
non-negative integers Ai > A2 > A3 > . . . such that ^ A^ = m, and likewise, let 

k 

/i = (/ii, yU25 /^3, • • •) be a partition of p. For fixed A,/i, the determinants of interest are 
given by 

^N-iif) = det{ax^-f,(-k+e)i<k,i<N- (A.l) 

For general A, fi, the determinant D^'^^{f) is not the minor of a larger Toeplitz matrix Tjy,j(/), 
but a minor with certain row and column indices shifted. For example if A = (2, 1,0,.. .) and 
/i = (1, 0, . . .), D^'^^i^f) is the determinant of from which we cross the third row and 

the second column, and then shift by 1 the indices of the Fourier coefficients lying on the 
first row. When Ai = fii, D^^_^{f) can be seen as a minor of the larger matrix Tat-i+AiI/) 
from which Ai rows and /ii columns are crossed out. 

The asymptotic values of the ratios D^^_^{f)/D]^_^{f) for large but fixed A and /i have 
been computed by Bump and Diaconis in [7]. Their expressions look complicated but are 
completely explicit and simple enough for small values of m and p. For two permutations 
TT and p, the numbers of fc-cycles in vr and in p are denoted respectively 7^ and 5k- If 
log = Yli Cne'""^, one defines 

n 

, ^ /fc^'=cr''=4!4?-''=^(-fcc.c_.) if7.>4, 
A(/,7r,p) = n I (A.2) 

k=i \k'-c'X'"'lk\Ll'~^''\-kcuc.u) if4>7fc, 



in terms of the Laguerre polynomials 

k=0 

Then one of the main results of [7i| asserts that, under suitable conditions on /, 



I — n \ / 



(A.4) 



where x'^(^) is the irreducible character of Sm associated to the partition A, evaluated at 
the group element tt, and similarly for x^{p)- 
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The previous result was proved in [3] under the conditions that the Fourier coefficients 
of log satisfy ^ |c„| < +00 and ^ |?t,| ■ |c„p < +00. For the apphcations we want to 

n n 

make of this result, the function / should be either /+ or given in (13.81) and (13.91) . and 
none of them satisfies these conditions, because the factor (2 — 2cos<y9)^^/^ has a root or a 
singularity on the unit circle. However the result (]A.4I) . written in the form given above, 
should hold even if these conditions are not satisfied since the right-hand side only depends 
on a finite number of Fourier coefficients of log /. 

Alternatively, and since the ratio is well-defined independently of the above two conditions 
on log/, one can circumvent the difficulty by regularizing the symbols /_|_ and /_, by defining 

^^1^) . (A.5) 
2t-2cos(pJ ^ ^ 

For t > 1, the two conditions are satisfied so that the formula flA.4|) can be used. Since 
the right-hand side is continuous at t = 1+, one can take the limit t ^ 1 at the end of the 
calculation, or equivalently use the Fourier coefficients Ck for the unregularized symbols /±. 

Let us now make contact with the minors needed in Section |Vl The ffist ones are given in 
(15. 9p . namely Mi(Ar;7v)(5' ± S')n+i/ det(S' ± S')n- In the present notation, they are equal to 
the minor Mi(jv.jv)(T/v(/)) divided by Dj^_i{f ) for / = f±. The matrix T/v(/) being Toeplitz 
of order + 1, its (A^, A^)-minor is equal to its (2, 2)-minor, itself equal to -D^^^]^(/) for the 
partitions X = fi = {1) of m = p = 1. The formula (]A.4[) yields 

Mi(^^^.)(T^(/)) Mi(2;2)(T^(/)) D'^'.,if) ^(0). . 
^™ F; T7\ = 1™ F; T7\ — = 1™ ~R T7\ = (-CiC_i) = 1 + CiC_i. 

(A.6) 

Taking f = f± and using (log/+)-|-i = V2 — 1 = — (log/_)±i yields the result given in (15.91) . 

The other minors we needed to compute are Mi(Ar;iv)(5' ± S")Ar+2/det(S' ± S")Ar+i in 
(I5.16p . The denominator is D^i^f) and the numerator, equal to Mi(3.3)(5' ± S')n+2 is noth- 
ing but D^'^(/), again for / = f±, and for the partitions A = yU = (1, 1) of 2. These two 
partitions label the alternating representation of 5*2, so that the above formula yields 

^ ^4°VciC-i) + jLf)(-2c2C_2)-J(c?c_2V_,C2)4^VciC-i)L«(-2c2C_2). 

(A.7) 

Taking f = f± and using the coefficients (log/+)±i = — (log/_)±i given above as well as 
(log/+)±2 = 3^2 - 4 = -(log/_)±2 yields the first part of flSTTj) . 
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The last minors to be evaluated in (15.161) are Mi(7v,Ar+i;Ar,Ar+i)(5' ± 5")iv+2/ det(5' ± S')]\r. 
It is not difficult to see that these ratios are equal to D^^^{f) / D j^_^{f) for / = /± and for 
A = yU = (2). Since these partitions label the trivial representation of 5*2, the formula (1A.4I1 
shows that this ratio is equal to the previous one in which all terms are taken positively, 
leading to 

lim -— = 

^Lf (-cic_i) + hf\-2c2C^2) + \{clc.2 + c2_iC2)4'Vcic_i)L«(-2c2C_2).(A.8) 
The values of the Fourier coefficients reproduce the second part of (15. 171) . 
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